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1 Introduction 



Nano-mechanical resonators can now be fabricated with frequencies approaching a 
gigahertz [1] . In this regime it is possible to consider using microwave transducers 
for their motion. Superconducting coplanar microwave cavities [2] offer a fabrica- 
tion geometry compatible with that used for nano-mechanical resonators and have 
recently enabled a new class of experiments in circuit quantum electrodynamics in 
the strong coupling regime [3J. In this paper we show how coherent states of a nano- 
electromechanical system (NEMS) may become entangled with the coherent states 
of a microwave field. The scheme may be used to improve the sensitivity for weak 
force measurement using nano-mechanical resonators. 

Oscillator coherent states are the quantum states that exhibit dynamics that is closet 
to the classical description of an oscillator. The correspondence extends to the case 
of two quadratically coupled oscillators initially prepared in coherent states as in 
that case the total system state, at all times factorises, into coherent states for each. 
Entangled coherent states however are highly non classical and can only arise when 
the Hamiltonian is higher order than quadratic in the canonical variablesjl]. Such 
states can be used to get an quantum enhancement in the signal-to-noise ratio for 
weak force detection^ or in optical quantum computing[H]. There are proposals 
for generating these optically [H |S]. Armour et al. [H] have suggested a scheme for 
putting a single nano-mechanical resonator into a superposition of coherent states. 

There are a number of ways that a nano-mechanical system may be coupled to a 
microwave cavity field. A possible experimental implementation is shown in figure 
[TJ The scheme we use here is based on that of Sun et al. PH] in which the interaction 
is mediated by an off-resonance interaction with a Cooper pair box (CPB) qubit. It 
is similar to the scheme presented by Zhang et al. [TT] . 

The coupling between a Cooper pair box charge system and the microwave field of 
circuit QED is given by[T2] 

H = AE C J2(N ~ n g (t)) 2 \N)(N\ - ^ £ (\N)(N + 1| + \N + 1)(N\) (1) 

N N=0 
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Cooper pair box 



Figure 1 : A nano- mechanical resonator is coupled to a microwave cavity via a Cooper 
pair box electric dipole. The nano-mechanical resonator carries a metal wire, is 
doubly clamped and suspended above a cut-away in the substrate. As it oscillates 
the DC bias voltage on a metal gate induces an AC bias voltage on the CPB. 

where 



with Cs the capacitance between the island and the rest of the circuit, C g is the 
capacitance between the CPB island and the bias gate for the island, and V g (t) is 
the total voltage applied to the island by the bias gate composed of a DC field, 
and microwave field in the cavity, v(t). Thus we can write V g (t) = Vg +v(t), where 
the hat indicates a quantisation of the cavity field. The Hamiltonian in Eq.([T|) is 
written in the Cooper pair number basis. In this basis the electrostatic energy of the 
first term is quite clear. The Josephson energy term describes tunneling of single 
Cooper pairs across the junction. This term is more traditionally (i.e. in mean-field 
theory) written in the phase representation as Ej cos 6. The connection between 
these two representations is discussed in [13] . 

We now allow for the possibility that the CPB (or its bias gate) is mounted on a nano- 
mechanical oscillator. This leads to a periodic modulation of the gate capacitance 



of the form 

C 9 (t) = C g {\ - x{t)/d) (2) 

where x is the displacement operator for the nano-mechanical oscillator and d is 
a typical length scale. Including both the time dependent cavity field and the 
modulated gate capacitance we can write 

n g {t) = nf + 5h g {t) (3) 

where 

r n {0) r 

sn 9 (t) = f e m - -j-m - ^fmt) (4) 

Using the usual restriction of the CPB Hilbert space to N = 0, 1, we can write the 
Hamiltonian as 

H = -2E c (l - 2nf)a z - ^-a x - 4E c 6h g (t)(l - 2nf - a z ) (5) 
where a z = |0)(0 - |1)(1|, a x = |1)(0| + |0)(1|. Define the bare CPB Hamiltonian 



as 



H CPB = -2E C (1 - 2nf)a z - ^-a x (6) 
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and diagonalise it as 



H C pb = 7} a z ( 7 ) 



where 



e=^Ej+[4E c (l-2N^W (8) 
and now the Hamiltonian takes the form, 

H = hjjj r a ] a + -a z - 4E c 5n g (t) [1 - 2n^ ] - cos 6a z + sin 6a x ] (9) 
2 y 

where we have now included the free Hamiltonian for the microwave cavity field, 

and 

9 = arctan [ y-r- ] (10) 

\AE C (1 - 2n?) J 

Operating at the charge degeneracy point, n g = 1/2 so that 6 = tt/2, the Hamil- 
tonian becomes 

H = ^ata + Wfe+l^-^r^^-^^-^x^^J^ (11) 
= fojj r a)a + hub^b + -a z — hg(a + (J)a x + h\(b + b^)a x + hx(a + a)){b + b\}&) 



where we have now included the free Hamiltonian for the NEMS (hvtfb), and 



hX = 2E C ^ (14) 
a 

= (15) 



where x rms = -JnJ (2mv) is the rms position fluctuations in the oscillator ground 
state. 

We now assume that the circuit resonance is at u T ~ 1 GHz, and that the circuit 
and the NEMS can be resonant. However, we detune the qubit from this resonant 
frequency by a few MHz [13]. We will make the rotating wave approximation and 
take the Hamiltonian in the interaction picture to be 

Hi = H6 a cra + hS(,¥b + hg(aa + + aV_) + hX(ba + + (16) 

with 5 a = u r — ujq and 5b = v — u are the detuning between the cavity resonance 
and the CPB and the NEMS and the CPB respectively, and hwo = e. 

We now assume that S a = 5^ = 5 (that is, u r = v) and that the time scale of interest 
is such that t » The effective interaction Hamiltonian is then 

H Ieff = H=[hx a ] a + Ml b ] b + hn (ab ] + a ] b) ] a z (17) 

with x = f2 = ^ and k = If the CPB is in a superposition of its two states, 
we can use this term to entangle coherent states of the NEMS with the cavity field 
as we now show. 



2 Non-dissipative case 

We first consider the case in which the NEMS, CPB and cavity are all non-dissipative. 
It will be shown that, in this case, unitary evolution under (I17jl results in entan- 
glement between the vibrational modes of the NEMS and the cavity field. We will 
assume that the NEMS and cavity field are in coherent states, |ao) a |A))f>) at t = 0. 
In the case that the CPB is in an eigenstate of a z , it can be shown that they remain 
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in coherent states at all times under (117j) and no entanglement can arise. However, 
if the CPB is prepared in an eigenstate of a x , entanglement is indeed possible. We 
first define cr z |±) = ±|±) and assume that the CPB is initially in its ground state, 
|— ), to write 

|tt(0)> =\a ) a \p ) b \-), (18) 
Enacting a n/2 pulse on the CPB puts it into the + eigenstate of a x , 

l*(0)> =^K)a|A)) 6 (|+) + |->). (19) 
The state at later times may be written 

W)> = 75 (|*+(0>«*l+> + l*-(*)>«*|-». (20) 

where 

\*±(t)) ab = e :fUH ±\a ) a ®\(3 ) b (21) 

with H = x a ^ a + Qtfb + n(a$ + a^b). Clearly the dynamics for a z = —1 is the time 
reversed dynamics for a z = +1. If however the CPB is in a superposition of the two 
eigenstates, we have the fascinating result of a superposition of the dynamics and 
the time reversed dynamics. 

Using the well known fact that the Hamiltonians H± transform coherent states into 
coherent states, the overall system state at time t can be written as, 

W)> = 75 (Ht))Mt)) b \+) + H-t))M-t)) b \-)) (22) 

where 

'k/Hq a (A-R)\ _ tRt/2 fa (A + R) k(3 



a(t) = e~ iuJt 
P(t) = e~ tuJt 



e iRt/2 + I "01"^^ _ \ e iR t /2 



iRt/2 



R 2R J \ 2R R 

koo /3 (A + R) \ _ im/2 _ f & (A - R) koq . 
R 2R r \ 2R + R 



(23) 
(24) 



and uj = (J] + x)/2, A = Q - x, R = VA 2 + 4k 2 . Note that for either 
case, a z = ±1, we have the conservation law |a(t)| 2 + |/3(t)| 2 = |a | 2 + |/?o| 2 - The 
vibrational modes of the NEMS and the cavity field have become entangled. To 
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detect the coherence implicit in this entanglement we need to readout a x on the 
CPB. 

The usual way to do this is to again subject the CPB to a 7r/2 rotation pulse, 
followed by a readout of a z for the CPB. This is equivalent to a direct readout of 
o x on the state in Eq. (j2*3l) . After the 7r/2 pulse, the system is then in the state, 

W{t)) = ~(Ht))Mt)) b -H-t))M-t)) b )\+) (25) 
+\ (Ht))Mt)) b + H-t))M-t)) b ) (26) 

It is then clear that a readout of the CPB state will project the cavity and NEMS 
into entangled coherent states, or 'cat states'. This is very similar to the NIST ion 
trap experiment to produce a cat state for the vibrational motion of the ion[16]. The 
coherence implicit in the cat appears in the probability for the CPB measurement 
results. For example, the probability to find the CPB in the state |— ) is 

P-(t) = ~ [1 + M ((a(t)\a(-t))mm-t)})] , (27) 



which may be written as 

-~ (Ht) - a(-t)\ 2 + W) - /3(-t)| 2 )] cos[$(t]} (28) 

where the interference frine; function of time are determined by the phase 

= 3{ a *(t)a(-t) + p*(t)P(-t)} (29) 

while the visibility of the interference is suppressed exponentially in the square 
of the distance between the coherent amplitudes and the time-reversed coherent 
amplitudes. This has a simple geometric interpretation given in figure El 

As a specific example consider the symmetric case in which Q = x — K so that 
A = 0. We further suppose that the nano-mechanical resonator is prepared with 
coherent amplitude j3 = B, taken as real, while the cavity is cooled to the ground 
state «o = 0. In this case we find that 

a(t) = -iBe~ iKt sin «t (30) 
P(t) = Be~ iKt cos Kt (31) 



so 




t)/2}} 



(32) 



In figure [5] we plot this function versus time for the case B = 4.0, k = 1.0. 
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Figure 2: A plot of the probability to find the CPB in the initial state |— ) as a 
function of t, with B = 4.0, k = 1.0. 

The collapse and revival of the interference can be interpreted geometrically, using 
a phase space representation, as shown in figure El Here we plot the orbits in the 
complex plane, parameterised by time, of the complex amplitudes a(±t) (on the 
left), and /3(±t) (on the right). The solid and dotted circles represent the minimum 
uncertainty circles for the coherent state. The dotted circle indicates a time at which 
we have complete overlap of the amplitude with its time-reversed form. At this point 
we expect from previous studies to see maximum interference |17j. This agrees with 
the result in figure [2j However, when the two coherent amplitudes do not overlap, 
the interference is exponentially suppressed and P-(t) is very close to 0.5. This 
geometric interpretation in terms of interference in phase space is analogous to the 
picture of Knight and co workers [18] for explaining the collapse and revivals in the 
Jaynes-Cummings model. 
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Figure 3: A graphical phase space representation of the generation of entanglement 
between the resonators and the CPB, and its manifestation as interference in P-(t). 
On the left (right) we represent the forward and reversed dynamics of the complex 
amplitude of a(6)-resonator. The parameters are the same as in figure El 

3 Dissipative case 

We now consider the case in which both the NEMS and cavity lose energy to 
their surroundings, which we treat as a zero temperature heat bath. This may 
be a reasonable assumption for the microwave cavity but is doubtful for a nano- 
electromechanical oscillator as they are currently operated, however it will at least 
capture the general features of decoherence as a best-case situation. We treat the 
damping in each case using the quantum optics master equation (weak coupling, 
rotating wave approximation). In the case of both oscillators approaching GHz 
resonant frequencies this approach is likely to be a reasonable approximation. In 
practice the decay rate for the nano-mechanical oscillator is likely to be greater 
than that for the microwave field. For simplicity we will ignore the qubit decay. 
The qubit itself is subject to both dephasing and spontaneous emission, both of 
which occur at rates of the order of 1 MHz[19]. As we now show, the rate of decay 
of coherence between the superposed coherent states is proportional to the modulus 
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square of the initial coherent amplitudes and thus likley to be very much greater 
than decoherence due to qubit decay. However, if necessary, the qubit decoherence 
can be included trivially as it simply causes an exponential decay, at fixed rate, of 
the off-diagonal term in the qubit basis. In general, the state of the system at any 
time is now described by a density matrix, 



Pit) = \ ( P++ (t)\+) (+| + p-(t)\-) (-1 + p + -(t)\+) (-1 + P- + (t)\-) (+|) (33) 
and the dynamics is given by the master equation, 

^ =-t[H I ,p}+ la V[a]p + lb V[b}p (34) 
where, for some operator A, 

V[A]p = ApA ] - \{A^ Ap + pA^ A) (35) 

and 7 a and % are respectively the dissipation rate constants for the cavity and the 
NEMS. 

In the case that the qubit is in an eigenstate of a z at t = 0, the resulting irreversible 
dynamics takes coherent states to coherent states, as does the reversible dynamics. 
This is a direct result of assuming zero temperature heat baths, and the state remains 
pure under irreversible evolution. We are thus led to adopt the ansatz that, as 
in the non-dissipative case, the cavity and NEMS will remain in coherent states, 
l a ±(^))a|/3±(^))fe|±) throughout unitary evolution. However, it is no longer true that 
the <j z = — 1 case is the time reversed version of a z = +1 case. We can then write 
for the diagonal terms, p ++ (t)\+) (+|, p (t)\— ) (— |, 

P++ (t) = \a + (t),f3 + (t)){a + (t),(3 + (t)\ (36) 
p __(t) = \a.(t),P_(t))(a.(t),p.(t)\ (37) 

where we have defined \ol,(3) = \a) a <g> \(5)b- 

Inserting this into (|3~41) and equating like terms gives us a system of ordinary differ- 
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ential equations for the coherent amplitudes: 



dt yi 

with initial conditions 



a± = =R(X"± + - ya± (38) 
/5± = T i(«a± + fi/3±)-y/3± (39) 
"« ± r+|/3 ± | 2 ) = -(7a|« ± r + 76 |/5 ± | 2 ) (40) 



a±(0) =a (41) 
/?±(0) = #,. (42) 

Note that a_(t) = a_(t)*, /?_(£) = /3+(t)*. 

We now turn to the off-diagonal term, p ^(t). Let us adopt the ansatz that p_ + has 

the form 

p_ + (t) = f(t)\a4t),^(t))(f3 + (t),a + (t)\ (43) 

with f(t — 0) = 1. The function f(t) controls the degree of coherence in the state 
due to the damping. For example, it is easy to see that the CPB coherence is given 
by 

<*«)(*) = »[/(*)] (44) 
In the protocol discussed for the non-dissipative case, it will determine P-(t) as 

P-(t) = \ [1 + » (f(t){a + (t)\a_(t)){P + (t)\P_(t)))} (45) 

By substituting Eq. (j43p into the master equation and equating like terms we arrive 
at the same system of differential equations as previously and in addition, 

ft ~ l f 7t ila - 12 + m " + |a+|2 + lP+]2) = f ^ aCe - a * + + lbf3 -^ (46) 

Using Eqs. (38-42) this may be written as 

I =/ew (47) 

where 

G(t) = 7 a (a-< - |«_| 2 ) + % {(3.(31 ~ I/ 9 - 1 2 ) • (48) 
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Using Eqs. (38-40) we can write 



Q{t) = 7 a(«4^)) 2 + l^-{t)f + j t (M*)| 2 + 1/9-WI 2 ) (49) 

So, 

f{t) = exp ((M*)| 2 + 1/3-WI 2 ) - (Kl 2 + |/3 | 2 ) + Jjt' la {a-{t')f + 7 6 (/3-(*')) 2 

(50) 

As the equations of motion for a±, (3± are a simple linear system, this may be easily 
computed by direct integration. The solutions are 

a_(*) = tf e -(7+-*("-t»/2))t + v e -^+- t{uJ+w/2))t (51) 

/?-(£) = Xe~ ( - 7+ ~ fi<w ~ u '/ 2 )' )< + ye~ ( - 7+_i ^ +u '^ 2 ^* (52) 

where 

[/ = ^[( M + A)/2 + i( 76 - To )/4]-^ (53) 

V = ^ [(to - A)/2 - z(76 - 7»)/4] + «^ (54) 

X = — [(to — A)/2 — z(7;, — 7a)/4] — k— (55) 

y = ^[(«; + A)/2 + i(76-7a)/4] + K^ (56) 

with w = ^Ak 2 - (j a - 7 6 + 2iA) 2 /4 and 7+ = (7„ + 7 6 ) /4. We then find that 

G{t) = la [x + {t)U 2 + x_{t)V 2 + 2UVy{t)] (57) 
+ 7b [x+(t)X 2 + x_(t)y 2 + 2Xyj/(t)] (58) 
+ (\a_(t)\ 2 -\a \ 2 ) + (\P-(t)\ 2 -\Po\ 2 ) (59) 

where 

J _ e -2t(7 + -i(c;=Fw/2) 

X± = 2( 7+ -i( WT »/2)) (60) 

1 _ g-2t(7 + -ioj) 

i/W = -57 n— (6i) 

2( 7+ - zu;) 

To gain an understanding of how decoherence works at short times, we consider 
the simple symmetric case: x — ^ = K ; la = lb = 7 with the initial conditions, 
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«o = 0, Po = B. In this case we find that 

a-{t) = ie'^-^B sin Kt 
P-(t) = e-^- lK)t B cos Kt. 



Thus 



f(t) = exp 



-(i 

2 v 



'jB 2 ( 1 — e ~7*+4iKt 



2 V 7 - 4m 
For short times (7^ << 1), this can be approximated by 

AB 2 



f(t) = exp 



7/A 3 + iB*K,yt 



2, -,,+2 



Thus the visibility of the interference decays at short times via the prefactor 

8B 2 



l/WI =exp 



-7« 2 t 3 



(62) 
(63) 

(64) 
(65) 
(66) 



a plot of which is given in fig. HI The cubic dependence on t has been seen in other 

; 




Figure 4: Decoherence as seen in the decay of \ f{t)\ 2 using the short time approxi- 
mation, with a = P, yt << 1, 7 = 0.1, x = f2 = 0, k = 1. 

models of this kind 1 17J. It is not unexpected. Coherence cannot be corrupted until 
it is created. Under unitary evolution, short time dynamics is always quadratic to 
lowest order in t, which then decays at a rate proportional to t for short times. The 
dependence of the decay rate on B 2 is also expected, as previous studies have shown 
that the coherence decays at a rate proportional to the square of the separation 
of the amplitudes of the superposed coherent states[T7j. Figl5] shows the general 
numerical solution for this case, where ao is again assumed to be real. Comparing 
figures (jU IHJ) we see that for these parameters the short time approximation is quite 
good. 
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Figure 5: Decoherence as seen in the decay of |/(£)| 2 using the full solution, with 
a = (3,yt « 1,7 = 0.1, x — Q — 0,k—1. 

4 Discussion and conclusion. 

We now make some optimistic estimates of the parameters that might be achieved in 
an experiment. From Walraff et al.[Hj we find that g ~ 6MHz and Ec/H ~ 5GHz. 
We will take the detuning between the two resonators and the CPB to be of the 
order of one MHz. We further suppose v = 1 GHz with a mass of m = 10 _21 kg [15] 
and thus x rms ~ 10~ 2 nm. With advanced fabrication it should be possible to achieve 
d ~ 20nm and thus A = 2(E c /h)(x rms /d) ~ 3MHz. It thus seems reasonable to 
expect that X/g = 0.5. We scale time so that x = 1, so that in dimensionless units, 
k/x — 0.5, Q/x — 0.25. The quality factor for the cavity resonance can be quite 
high, perhaps as high as 10 4 . With a resonance frequency of 1GHz, this gives in 
dimensionless units that 7 = 10~ 3 . The quality factor for the mechanical resonance 
is not as good, so we take it to be a factor of 10 smaller so that 7^ = 10~ 2 . This 
might be achievable for a very small nano-mechanical resonator as considered here. 
Let us take the microwave resonator and mechanical resonator to have low coherent 
amplitudes Oq = /3q = 2.0, which should be achievable at milli Kelvin temperatures. 
In figure ([U]) we plot the probability to find the CPB in the ground state using the 
protocol of section [2j Interference effects are clearly evident giving some hope that 
the quantum entanglement implicit in this result could be observed. 

We have presented a method to entangle the coherent motion of a nano-mechanical 
resonator with a coherent state of a microwave cavity field. In so far as the nano- 
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Figure 6: The probability to find the CPB in the ground state using the protocol 
in section [2] including dissipation, with a = (3 — 2, 7 a = 0.001, 7& = 0.01, x — 
1.0, Q = 0.25, k = 0.5. 

mechanical resonator is a mesoscopic mechanical system, this is already a fascinating 
prospect similar in many ways to the two-slit interference of a fullerne molecule [20j . 
However, it may also be of some use. Entangled coherent states can be used to 
improve the sensitivity for weak force measurements [5]. As weak force detection 
is one of the primary motivations for building nano-mechanical systems, this may 
prove a useful way to improve performance. 
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Figure captions 



1. A nano- mechanical resonator is coupled to a microwave cavity via a Cooper 
pair box electric dipole. The nano-mechanical resonator carries a metal wire, 
is doubly clamped and suspended above a cut-away in the substrate. As it 
oscillates the DC bias voltage on a metal gate induces an AC bias voltage on 
the CPB. 

2. A plot of the probability to find the CPB in the initial state |— ) as a function 
of t, with B = 4.0, k = 1.0. 

3. A graphical phase space representation of the generation of entanglement be- 
tween the resonators and the CPB, and its manifestation as interference in 
P-{t). On the left (right) we represent the forward and reversed dynamics of 
the complex amplitude of a(6)-resonator. The parameters are the same as in 
figure El 

4. Decoherence as seen in the decay of \ f{t)\ 2 using the short time approximation, 
with a = (3, jt « 1,7 = 0.1, x = fi = 0, k = 1. 

5. Decoherence as seen in the decay of \f(t)\ 2 using the full solution, with a = 
/3,7f«l,7 = 0.1, X = n = 0,K=l. 

6. The probability to find the CPB in the ground state using the protocol in 
section [5] including dissipation, with a — j3 — 2, 7 a = 0.001, 75 = 0.01, x — 
1.0, Vl = 0.25, K = 0.5. 
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